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Abstract 

The logical structure and the basic theorems of time-dependent current density functional the- 
ory (TDCDFT) are analyzed and reconsidered from the point of view of recently proposed time- 
dependent deformation functional theory (TDDefFT). It is shown that the formalism of TDDefFT 
allows to avoid a traditional external potential-to-density/current mapping. Instead the theory is 
formulated in a form similar to the constrained search procedure in the ground state DFT. Within 
this formulation of TDCDFT all basic functionals appear from the solution of a constrained univer- 
sal many-body problem in a comoving reference frame, which is equivalent to finding a conditional 
extremum of a certain universal action functional. As a result the physical origin of the univer- 
sal functionals entering the theory, as well as their proper causal structure becomes obvious. In 
particular, this leaves no room for any doubt concerning predictive power of the theory. 
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I. INTRODUCTION 



Last decades witness a growing popularity of the time-dependent density functional the- 
ory (TDDFT) l|] as a practical tool for studying dynamics of various quantum many-body 
systems (see {2I and references therein). The range of applications of TDDFT (including 
its current-based version, the time- dependent current density functional theory (TDCDFT), 
3|,|j]) is impressively broad. Interaction of atoms, molecules and solids with electro- magnetic 
radiation, transport phenomena in nanostructured systems, dynamics of nuclear matter and 
ultracold trapped atomic gases form by far not a complete list of modern TDDFT/TDCDFT 
applications. 

The first fundamental contribution to TDDFT, which can be also considered as a birth 
of the theory, was made by Runge and Gross in 1984 Many importand developments 
and generalizations of TDDFT which appeared in the last 25 years deepened significantly 
our understanding of this approach (a good discussion of the current state of a formal 
theory can be found in a review [5| and in Chapters 1-7 of the book [2]). Nonetheless 
surprising observations, fundamental puzzles as well as misunderstandings related to formal 
aspects of TDDFT regularly appear up to now. The most known and, possibly, historically 
the first puzzle of this kind was a so called causality paradox closely related to problems 
with variational definition of xc potentials in TDDFT. Presently we have many various 
explanations and resolutions of this paradox , but, surprisingly, the most 

straightforward and elementary resolution of the problem appeared only in the last year 
I Anothj puzzling =u.pH=e that appealed m the last yea. was and obse™tio„ n,ade 
by Baer [10|. He has found that in a lattice formulation of the theory, in contrasts to 
the ground state DFT, there are apparently reasonable time-dependent densities which are 
not v-representable. This problem was discussed and clarified in two very resent papers 
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12|, but it became quite clear that v-representability, which is commonly taken for 



granted in TDDFT and TDCDFT, can be a serious issue. The second fundamental issue of 
TDDFT, the density-to-potential mapping, was also under debate in the last year 
The authors of Ref. [l^ attempted to analyze the causal structure of the mapping from 
the density to external potential and arrived to the conclusion that TDDFT in the Kohn- 
Sham (KS) formulation can not predict the dynamics of the density as the xc potential 
does depend on a "future" in a form of a second time derivative of the density. In a 
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recent comment Maitra, et. al. '13!] have demonstrated that this conclusion is based on a 
misunderstanding of the basic mapping theorem 



15|, and proved that KS-TDDFT has a 
full predictive power. However at the present stage of the theory this can be shown only 
for the KS-TDDFT, while for the other, for example, hydrodynamic formulations [la] 
of the theory the question raised in Ref. is still formally valid. In general we should 
honestly say that the causal structure of all basic functionals in TDDFT and TDCDFT is 
not well studied and understood. Hence in the last year we have seen the whole bunch of 
papers debating actually all most fundamental points of TDDFT and related approaches. 
Therefore it is timely to carefully reexamine the foundations and the logical structure of 
the theory and formulate it, from the very beginning, in a way that takes into account our 
current level of understanding and the most recent developments in the field. 

This work is an attempt of such reexamination. Below I concentrate of the current den- 
sity version of the theory, having in mind its generality and numerous recent applications of 
TDDFT-based methods to the transport theory where the exact knowledge of the current 
is vital. From conceptual point of view TDCDFT is a reduced theory which allows one 
to describe the behavior of the current j and the density n formally ignoring the rest of 
complicated microscopic dynamics of a quantum many-body system. In other words TD- 
CDFT is a closed theory of convective motion. Therefore it is most natural to formulate 
such a theory in a form of a closed system of equations of motion for the basic variables, 
the density and the current. This is the hydrodynamic formulation of TDCDFT, which 
analyzed in detail in Sec. [Tll I show how hydrodynamic equations of motion are related to 
the microscopic many-body dynamics, and how one can make them closed by employing 
the traditional formulations of the TDCDFT mapping theorem 3], Q]. The main concern 
of Sec. [TTl is to prove an equivalence of using different sets of basic variables which can be 
used in TDCDFT to describe the convective degrees of freedom. These are the current and 
the density, the velocity field and the density, (v, n), in the Eulerian description of 

convective dynamics, or a continuous set of trajectories x(^, t) of infinitesimal fluid element 
in the Lagrangian description. The hydrodynamic equations of TDCDFT for all three pos- 
sible choices of basic variables are derived and compared. 1 also derive a KS formulation 
of TDCDFT and show how the basic universal functionals entering the hydrodynamic TD- 
CDFT and KS-TDCDFT are related to each other. The most important part of the paper 
is Sec. mil where the foundations of TDCDFT are reconsidered from the point of view of 
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recently proposed time- dependent deformation functional theory (TDDefFT) 17|. I prove 
the equivalence of two theories and present a new logical structure of TDCDFT which is 
based on TDDefFT. The idea is very simple and physically transparent. From the very 
beginning we take the Lagrangian trajectories as a set of basic variables for a future closed 
theory of the convective motion, and separate the convective degrees of freedom by trans- 
forming the many-body theory to a reference frame moving with fluid elements. After this 
transformation an appearance of a closed theory of the convective motion becomes obvious. 
In fact, TDDefFT and thus TDCDFT take a 
search formulation of the ground state DFT 18 



brm which very similar to the constrained 



In this paper I am trying to avoid long algebraic manipulations and simply state the re- 



sults when they are physically plausible (some detailed calculations, especially 



those related 



formalities of nonlinear transformations of coordinates can be found in [16|, Il7|, |2l|. Instead, 
I am concentrating on the logical structure of the theory and on discussions of some delicate 
issues, such as equivalence of different descriptions, uniqueness and existence of solutions, 
etc. I hope the this way of presentation is more suited to the purpose of this work. 

II. HYDRODYNAMIC FORMULATION OF TDCDFT 

A. Starting point: General formulation of the quantum many-body problem 

Let us consider a system of identical particles in the presence of time dependent 
external scalar f/(x, t) and vector A(x, t) potentials. The corresponding many-body wave 
function \I'(xi, . . . ,XAr,t) is a solution to the time-dependent Schrodinger equation 



Z(9t^(Xi, . . . , X^v, t) = i^'^(Xi, . . . , Xtv, t) 



(2.1) 



with the following Hamiltonian 

N 



Z(9x. + A(x,-,t))2 



2m 



i5^V(|x,-x,|) 



where y(|x — x'|) is the interaction potential. For a given initial condition, 

^(Xi, . . . , Xat, 0) = ^o(xi, . . . , Xat), 



(2.2) 



(2.3) 



the dynamics of the system is completely specified by Eq. (12.11) 



In most physically important situations it is not necessarily to know the full many-body 
wave function. Normally the experimentally measurable response of the system to external 
probes can be described in terms reduced "collective" variables - the density of particles 
n{'x,t), and the density of current j(x, t) 

n(x,t) = p(x,x,t), (2.4) 

j(x,t) = — lim(9x-9xOp(x,x',t)--A(x,t), (2.5) 
zm x'— >x m 

where p(x, x', t) is the one particle reduced density matrix 

N 

p(x, x',t)=N Yl c^Xj^*(x, X2, . . . , x;v, t)^(x', X2, . . . , X^, t) (2.6) 

The main idea of TDCDFT is to reduce, at the formally exact level, the problem of 
calculation of the density and the current to solving a closed system of equations which 
involve only n(x, t) and j(x, t). In the next two subsections we describe a few equivalent 
ways to formulate such a closed theory. 



B. Local conservation laws and TDCDFT hydrodynamics in Eulerian formulation 

Using the microscopic definitions of Eqs. (12. 4p and (12.51) . and the Schrodinger equation 
(12. ip one can derive the following hydrodynamic equations of motion the density and the 
current 

dtn + d^^j^ = 0, (2.7) 
mdtj^ - [j X B]^ - nEf, + d^.U^,^ = (2.8) 

where E(x, t) and B(x, t) are electric and magnetic fields generated by the external time- 
dependent scalar and vector potentials 

E(x,t) = -9tA(x,t) -9xf/(x,t), (2.9) 
B(x,t) = d^x A(x,t). (2.10) 

Equation (12. 7p is the usual continuity equation, i. e., a local conservation law of the number 
of particles. The equation of motion for the current, Eq. (12.80 . physically corresponds to 
a local momentum conservation law (or a local force balance): the time derivative of the 
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current equals to a sum of the external and internal forces. Importantly, the local internal 
force (the last term in Eq. fl2.8p ) has a form of a divergence of a second rank tensor. Therefore 
the bulk internal force vanishes after a volume integration, as it is required by the Newton 

ow tensor 11^,^ entering Eq. fl2.8p contains a kinetic and an 



third law. The momentum f 
interaction contributions 21 



22| . U^^ = nj^™ + n™*, which are expressed in terms of the 



many-body wave function as follows 

"2 



hm iP;Pl + P:Pl)p{^. x', t) - ^dlni^, t) 



(2.11) 

nj*(x,t) = -i I rfx'^^^^^'rfAr(x + Ax',x-(l-A)x,t), (2.12) 



where = —idxi^ — A^(x, t) is the kinematic momentum operator, and r(x, x',t) is a two 
particle reduced density matrix 

N 

r(x, X', t) = NiN - 1) / n ^Xj^*(x, X', X3, . . . , X^, t)^(x, X', X3, . . . , XAT, t). (2.13) 

j=3 

A formal possibility to formulate a closed theory for calculation the density and current 
distributions follows from the mapping theorem of TDCDFT j^, J]. The main statement of 
this theorem can be formulated as follows. For a given initial state \l/o a map of the external 
potentials to the current, {U, A) H-i> j, is invertible and unique up to a gauge transformation, 
provided the potentials are analytic in time, and the current density is v-representable. 
(Note that here the term "v-representability" is understood in a broad sense: a current j is 
called v-representable if it can be produced by some external 4-potential {U, A)). 

In fact, the above mapping theorem states that the potentials are unique (modulo a 
gauge transformation) functionals of the initial state and the current density, U = f/[\I^O; j] 
and A = A[^Oi j]- This immediately implies that the many-body wave function \l/(t) and, 
therefore, any physical observable is also a functional of \I^o and j. In particular, inserting 
the functionals A[^o, j] and ^ = ^[^^o, j] into the definitions of Eqs. fl2?TT|) and (I2TT2D we get 
the exact momentum flow tensor as a unique functional of the initial wave function and the 
current II^i, = n^i,[\l/o, j]. This functional is universal in a sense that it does not explicitly 
contain the external potentials, but is uniquely recovered from a given current and an initial 
state (in the following for the sake brevity we omit \l/o in the arguments of the functionals). 
Substituting the functional n^j,[j] into Eq. (12.81) we obtain a closed system of equations of 
motion for n(x, t) and j(x, t). Hence from the system of Eqs. (I2.7l) - (l2.8p we can in principle 



determine the dynamics of the density of particles and the density of current, avoiding, at 
least formally, the explicit solution of the full many-body problem. TDCDFT represented in 
a form of the closed system fl2.7l) - (l2.8p can be viewed as an exact quantum hydrodynamics. 
It is worth noting that the hydrodynamic formulation of TDCDFT is analogous to the 
brmulation of the static DFT in a form of the direct Hohenberg-Kohn variational principle 



23|. 



A connection of the TDCDFT hydrodynamics to the standard mechanics fluids [24'] can 
be made more obvious if we switch the basic variable from the current j to the velocity field 
V = j/n. It is also useful to extract from the full momentum flow tensor II^i, its exactly 
known part - the flow of momentum due to convective motion of the fluid, mnv^Vy, 

li^^ = mnv^Vy + P^y (2.14) 

where P^^ is the stress tensor which is responsible for a local internal force related to a 
relative motion of particles "inside" an small moving fluid element (see Sec. IIIII below for 
a more detailed discussion). Using the representation (12.141) and expressing all currents in 
terms of the velocity field we transform equations of motion fl2.7p and (12.81) to the following 
standard "Navier-Stokes" form 

{dt + v9x)n + nd^p^v^ = 0, (2.15) 
midt + vd^)v^ - [v X B]^ -E^ + ^d^.P^,[v] = 0. (2.16) 

Since the map (n,j) i-^ (^, v) is one-to-one we are allowed to replace the functional depen- 
dence of the stress tensor on the current by the functional dependence dependence on the 
velocity field, P^u[ji\ ^ Pi^iul^]- It is also worth noting the knowledge of only the current j 
or the velocity v is sufficient to recover the density by integrating the continuity equation 
of the form of (12. 7p or (12. 7p , respectively. 

C. Kohn-Sham construction in TDCDFT 

Practical applications of any DFT always rely on the Kohn-Sham (KS) construction 
25| . which in the present time-dependent setting can be introduced as follows (see, e. g., 
161]). Let us consider a fictitious system of noninteracting particles in the presence of an 
electromagnetic field generated by the external 4-potential {U, A), and by some selfconsistent 



vector A^*^ and scalar JJ = JJ^^ + U potentials, where U is the usual Hartree potential. 
The dynamics of this system is described by a set of one-particle Schrodinger equations for 
KS orbitals 0j(x, t), j = 1 . . . N 

tdt(p, = ^{td^ + A + A^'f<j)j + {U + f/^^«)</), (2.17) 

Obviously the density ns and the velocity V5 of the KS system satisfy the continuity equation 
( I2.15P and the force balance equation of the form of Eq. f l2.16p . but with the stress tensor 
P^u, and the external Lorentz force being replaced, respectively, by the kinetic stress tensor 
T^j, of the noninteracting KS particles, and by the Lorentz force corresponding to the total 
effective 4-potential. From the requirement that the KS density and current reproduce the 
density and the current in the real interacting system we get the following selfconsistency 
equation connecting the xc potential to the stress tensor functional 

dtA^;^' - [v X (^x X A^^)]^ + ^..f/^'^" = -i9,.AP^,[v], (2.18) 

where AP^j^[v] = -P;x;/[v] — T^j,[v] is a difference of the stress tensors in the interacting and 
noninteracting KS systems. Equation (12.181) can serve as a most general definition of the 
xc potentials. For a given stress tensor functional (the right hand side) it defines the xc 
4-potential {U^^, A^'^) up a gauge transformation. 

It is important to stress out that the KS construction is only an auxiliary formal device 
for solving the general collective variable theory in a form of closed equations of motion, 
Eqs. (12.71) and (12. 8p ( or, equivalently, Eqs. (12.151) and (I2.16P ). for basic variables - the 
current/velocity and the density. The situation is very similar to the static DFT where the 
KS costruction serves merely as useful mathematical trick for transforming the fundamental 
Hohenberg-Kohn variational principle to a system of differential equation for one-particle 
orbitals. 



D. Does TDCDFT/TDDFT have a predictive power? 



as TDDFT) looks 



14| . Let us assume 



Despite conceptually the traditional formulation of TDCDFT (as wel 
clean and very similar to the static DFT, it may cause some confusions 
that the stress tensor functional is known. The main goal of the theory is to calculate the 
current and the density at any instant t by propagating the equations of motion (I2.7p and 
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(12. 8p starting from the initial time, provided the initial conditions for j and n are given. 
Apparently the propagation of Eq. (12. 8p is only possible if n^,^[j] is a retarded functional 
of the current. In other words at a given instant t the stress tensor should only depend 
on the currents j(t') at previous times t' < t. In particular, it should not contain local 
in time terms which depend on the time derivatives (of any order r > 1) of the current. 
However, both the original Runge-Gross proof and its generalization to TDCDFT by 
Ghosh and Dhara j4| establish only the uniqueness of the maps j i-^ A and j i-^ \E' , but 
do not state anything about the properties of functionals A[j] and ^E'p]. On the other 
hand if we consider the local momentum conservation law, Eq. (12. Sp . and interpret it as 
an equation which determines the external 4-potential for a given current, we immediately 
observe that the external force, and thus the external potential considered as functional 
of the current does indeed contain a local in time term proportional to Taking this 
fact naively one may conclude that the functionals and n^,^[j] should also contain such 
unwanted terms. Therefore the propagation of the basic equation Eq. ( 12. Sp (or equivalently 
Eq. (I2.16P ) becomes problematic. Hence it appears that TDCDFT/TDDFT does not have 
a real predictive power. 

In a recent paper Maitra et. al. clearly demonstrated that the above conclusion is incorrect 
for the KS formulation of TDDFT |13]. In fact, this is obvious from the van Leeuwen proof 
of the TDDFT mapping theorem jisl, a nd its generalization to TDCDFT by Vignale j^. 



The construction proposed in Refs. 3 



Iq shows that the difference of potentials driving the 



dynamics two different systems with the same current (e. g., interacting and noninteracting 
systems in the KS case) is uniquely recovered from the given density /current and has a 
proper causal structure. Unfortunately this methodology does not allow, at least directly, 
to make any definite statement about the retardation properties of the many-body wave 
function as a functional of the current and, what is much more important, about the 
causal structure of the stress tensor functional entering the collective variable theory defined 
by Eqs. (I2.15P and (I2.16p . In Sec. Illll we show how these questions are resolved within a new 
constructive reformulation of the theory based on TDDefFT. As a first step in this direction 
we discuss one more alternative form of the TDCDFT hydrodynamics. 
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E. TDCDFT hydrodynamics in the Lagrangian form 



In general TDCDFT is a closed formalism that allows one to describe a convective motion 
of a quantum many-body system driven by an external field. Usually the convective motion 
is characterized by the density of particles n(x, t) and the density of current j(x, t) or a 
velocity field v(x, t). An alternative way to completely characterize the convective motion is 
commonly referred to as a Lagrangian description. Let us consider the system as a collection 
of infinitesimal fluid elements (so called "materials"). Every fluid element can be uniquely 
labeled by a continuous variable $, - its position at the initial time t = 0. The Lagrangian 
description can be viewed as tracking the motion of those infinitesimal elements of the fluid. 
In other words the convective motion of the system is characterized by a (continuous) set 
of trajectories x(^,t), where the argument $, indicates the starting point of the trajectory 
(the unique label of the element). Let us show that in the case of dynamics of a quantum 
many-body system the map v(x, t) x(^,t) is unique and invertible. 

For a given velocity v(x, t) the Lagrangian trajectory is solution to the following Cauchy 
problem 

dMt t) = v(x(^, t),t), x(^, 0) = t (2.19) 

It is known (see, for example, Ref. 26j) that Eq. fl2.19p has a unique solution x(^,t) if 
the function v(x, t) is continuous and satisfies the Lipschitz condition in spatial variables 
(i. e. there exists a constant L > 0, such that ||v(x) — v(x')|| < L ||x — x'|| for any x and 
x'). Apparently the velocity field coming from a physical wave function does satisfy these 
requirements. Physically the Lipschitz condition prevents generation of folds and shock 
fronts which are clearly absent at the microscopic level for the Schrodinger dynamics. Hence 
from a given physical velocity field we uniquely construct a set of Lagrangian trajectories. 
Every trajectory x(^, t) is uniquely determined by its initial point ^, which means that given 
the initial position £ of a fluid element we can always find its coordinate x = x(^, t) at any 
instant t, and, tracing the trajectory in a reverse order, from a given position x at time t one 
uniquely recovers the initial point of the trajectory, ^ = ^(x, t). In a more formal language, 
the map ^ i— x : x = x(^,)f:) is unique and invertible. This fact allows us to recover the 
Eulerian variables, i. e., the velocity field and the density, from given Lagrangian trajectories 



v(x,t) 



dt 



(2.20) 
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where ^(x, t) is the inverse of x(^,t), 'n,o(x) is the initial density, and ^/ g{^, t) = J(^,t) = 
det(^) is the Jacobian of the transformation of coordinates x — ^. Equation fl2.20p is an 
obvious consequence of Eq. (12.191) . while Eq. fl2.2ip is a direct solution of the continuity equa- 
tion (12.151) . Therefore the function x(^,t) indeed completely characterizes the convective 
motion of a system. 

The basic equation in the Lagrangian description of collective dynamics is the equation of 
motion for a fluid element. This equation can be straightforwardly derived from the equation 
of motion for Eulerian velocity v(x, t), Eq. (12.161) . by making a transformation coordinates 
X — s> ^, i. e., by considering the initial points ^ of Lagrangian trajectories as independent 
spatial coordinates. Under this transformation the convective derivative, dt + v9x becomes 
simply dt so that the first term in Eq. (I2.16P transforms to mx(^, t), while the divergence of 
the stress tensor in the last term in Eq. (I2.16P becomes a covariant divergence in the space 
with metrics g^u{$,,t) induced by the transformation from x- to ^-coordinates. Hence after 
the transformation of coordinates we arrive at the following equation of motion for a fiuid 
element 

mx'^ - E^(x, t) - [i X B(x, t)l + V,P^^[x(£, t)] = 0, (2.22) 

where P^^i^^i) is the original stress tensor Pa/}(x, t) transformed to the new coordinates 



according to the standard rules 



27| 



P,Att)=g^^Pa,{^{U),t) (2.23) 
The nabla-operator in Eq. (12.221) stands for a covariant divergence 

V.P; = ^%v^P; - (2.24) 
and the metric tensor in the ^-space of "initial positions" is defined as follows 

(2.25) 

The equation of motion (I2.22p has to be solved with the initial conditions x(^, 0) = ^ and 
^(^) 0) = '^o(€)) where vo(x) is the initial velocity distribution calculated from the initial 
many-body wave function of Eq. (12. 3p . 
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The first three terms in Eq. fl2.22p correspond to a classical Newton equation for a particle 
moving in the external electromagnetic field, while the last, stress term takes care of all 
complicated quantum and many-body effects. Because of the uniqueness and invertibility 
of the map v(x, t) x(^,t) the transformed stress tensor can be considered as a unique 
functional of the Lagrangian trajectories, P^i, = P^,^[x(^, t)]. Hence Eq. f l2.22p is a closed 
equation of motion, which, at the formally exact level, completely determines the collective 
dynamics of the system. This is the basic equation of TDCDFT in the Lagrangian form. 

From the first sight the representation of TDCDFT hydrodynamics in the Lagrangian 
form of Eq. (12.221) does not bring anything fundamentally new. This is indeed true if 
one follows a route outlined in this section: starting from the traditional formulation of 
the TDCDFT mapping theorem and via the Eulerian equation of motion for the density 
and the current. However, in the next section we will see that using the ideas of the 
Lagrangian description one can reformulate the whole theory, including all basic theorems, 
in a constructive way that also ends up with the equation of motion fl2.22p . but provides 
us with a clear constrained search-like procedure for calculating the basic stress tensor 
functional. 

III. TIME-DEPENDENT DEFORMATION FUNCTIONAL THEORY 
A. Many-body theory in a comoving reference frame 

Conceptually TDCDFT is a reduced theory aimed at describing only the convective 
motion of the system without a detailed knowledge of the full dynamics of all microscopic 
degrees of freedom. Therefore it looks natural to start the consruction of such a theory by 
separating the "convective" degrees of freedom at the very beginning, i. e. at the level of the 
full microscopic many-body theory. The Lagrangian description is perfectly suited for this 
purpose. Indeed, in this formalism the convective dynamics is characterized by the motion 
of fluid elements. Therefore it can be easily separated from the microscopic dynamics of 
quantum particles by transforming the many-body theory to a local noninertial reference 
frame moving along the Lagrangian trajectories. 

At the formal level one proceeds as follows. Consider a reference frame defined by some 
(unspecified for the moment) velocity field v(x, t) which is required to be continuous and 
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Lipschitz in spatial variables. By solving the Cauchy problem of Eq. (12.190 with the above 
velocity in the right hand side we get the local trajectories x(^,t) of the frame. The trans- 
formation of the theory to the new reference frame corresponds to a transformation of 
coordinates Xj ^j, with x(^,t) being the transformation function, i. e. Xj = x(^j,t), in 
the many-body Schrodinger equation (12. ip . It is convenient to define the many-body wave 
function \E'(^i, . . . , ^at, t) in the new frame as follows {l3| 

N 

*(^i, ■■■,^N,t) = l[g'* (^„ t)e-^^-(^-*)vl>(x(^i, t), . . . , x(^;v, t),t), (3.26) 

where >S'ci(^,t) is the classical action of a particle moving along the trajectory x(^,t) 

SciiC t) = 1^ [jim? + i(t) A(x(t), t) - t/(x(t), t)] . (3.27) 

Equation (I3.26P is a relatively straightforward generalization of the transformation to a 
homogeneously accelerated frame, which is used, for example, in the proofs of a harmonic 
potential theorem [28|, l29| . The exponential prefactor accounts for the phase acquired due 
to motion of the frame, while the factor Y[j'^=i9^i^j^^) aimed at preserving the standard 
normalization of the wave function (^l/|^E') = 1 after a non- volume-preserving transformation 
of coordinates. 

Performing a transformation of coordinates, Xj — ^j, in Eq. (12.10 . and using the defini- 
tion (I3.26P we obtain the many-body Schrodinger equation in the frame moving with some 
velocity v(x, t) 

tdt^i^u ...,^N,t) = H[g,,, A]^{iu ...,^N,t) (3.28) 
The Hamiltonian in the new frame takes the form 

N ^ . — - fiu ^ _ 

i=i Mi 

where Kj^^j, = —id^t^ — A^{$,j,t), g^^ = g^'^{^j,t), and l^^,^. is the distance between jth and 
kth particles in the moving frame (the length of geodesic connecting points $,j and $,k in the 
space with metric g^u)- An effective vector potential A.{^,t) entering the Hamiltonian is 
given by the following equation 

A, = —x^ + ^^.(x,t) - d,.S.,{tt). (3.30) 
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Physically this effective vector potential describes a combined action of the external and 
inertial forces in a local noninertial frame. Since we defined the moving frame in such a way 
that at if: = it coinsides with the laboratory one (this formally follows from the condition 
x(^) 0) = ^) the initial condition to the transformed Schrodinger equation fl3.28p remains 
unchanged: 

$(^i,...,^jv,0) = vl>o(^i,...,^^). (3.31) 

Equations fl2.22l) - fl3.31l) completely determine dynamics of the quantum iV-particle system 
in an arbitrary local noninertial frame. 

Since our aim was to reformulate the theory in a particular frame moving with a physical 
flow (this frame is called comoving, or Lagrangian) we need to impose an additional local 
condition to specify the required frame. By definition in the comoving frame the current 
density is zero everywhere and at all times, while the density of particles stays stationary 
and equal to the initial density distribution no(^). Hence the most natural frame- fixing 
condition is the requirement of zero transformed current density j(^,t) = 0. Explicitly this 
condition reads 

Ait t) = hm (9^ - d^,m, t) (3.32) 

where p(^, t) is the one particle reduced density matrix calculated from the transformed 
wave function 




(3.33) 



The frame-fixing condition (13.321) simply states that in the comoving frame a "paramagnetic" 
current (the right hand side of (13.321) ) is exactly cancelled by the "diamagnetic" contribution 
(the left hand side of ([332])). 

The Schrodinger equation (13.280 . the definition of the effective vector potential (13.301) . 
and the zero current condition (I3.32p constitute a closed system of equations that deter- 
mine dynamics of the many-body system in the comoving reference frame. In principle one 
can eliminate the effective vector potential from this system by substituting yA.{$,,t) from 
Eq. (I3.32P into Eqs. (I3.30p and (13.280 . As a result we will get a system of two first order (in 
time) differential equation for two functions - x(^, t) which describes the convective motion 
on the system, and \E'(^i, . . . , ^tv? ^) describing the rest of microscopic degrees of freedom in 
the frame moving with the convective flow. The equations have to be solved with the initial 
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conditions Eq. (I3.3ip for and x(^, 0) = ^ for the Lagrangian trajectory. The system of 
Eqs. f l3.28p . fl3.30p and (13.320 is equivalent to the original linear Schrodinger equation (12.11) . 
Therefore we can guarantee that there exists a unique solution of the corresponding initial 
value problem. 



B. A "constrained search" formulation of TDDefFT 



Let us discuss possible procedures for solving the system of Eqs. (13.280 . (13.300 and (13.320 . 
but, before doing that, it makes a sense to rewrite it in a more physical and clear form. 

First of all, we note that because Eq. (13.301) contains the classical action Sci{^,t) of 
Eq. (13.270 . it is nonlocal in time. It is convenient to remove this nonlocality by differentiating 
Eq. (I3.3O0 with respect to t. Taking the derivative we reduce Eq. (13.301) to the following 
form 



9e 



I If- 

f / . . I ^ 1,1 — I — lA. lai.^ I. II.. — I — , 



mx" = ^^(x, t) + [i X B(x, t% + TTTr^tA, (3.34) 



which is exactly the classical Newtonian equation for a particle moving in the external 
electromagnetic field and, in addition, in the "electric" field, generated by the effective 
vector potential. 

Comparing Eq. (I3.34p with Eq. (I2.22p we immediately observe that these two equation 
would be identical if the time derivative of the effective vector potential would equal to the 
covariant divergence of the stress tensor. To see that this is indeed the case we consider the 
local force balance equation that follows from the many-body Schrodinger equation (I3.28p 
in the local noninertial frame. Apparently this equation is of the form of Eq. (12.80 . but with 
the usual divergence of the momentum flow tensor replaced by covariant divergence of the 
stress tensor in the ^-space. Namely, 

^Jm - ?(% A - %^,) + ndtA^ + V^V.P; = 0. (3.35) 

where the stress tensor in the s pac e with metric g^y, can be conveniently determined from 
the following universal formula 17 , [21I, 3^ 

h^{i.t) = - {^\K[9a,,A]\^) (3.36) 

with the Hamiltonian defined by Eq. (I3.29p . An explicit form of the stress tensor operator 
P^y entering Eq. (13.361) can be found, for example, in Ref. 2l|. At the moment the important 
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point is only that P^viQap-, ^] is an explicitly known and local in time functional of the metric 
tensor and the effective vector potential. Since in the comoving frame the transformed 
current density j is zero, only the last two terms survive in Eq. (13.351) . Therefore in our 
frame of interest the force balance equation reduces to the following identity 

dtA, = -^V,P; (3.37) 

Inserting this identity into Eq. (13.341) we exactly recover the basic equation of hydrodynamics 
in the Lagrangian description, Eq. (I2.22p . The important progress is that the stress tensor 
in this equation is now defined entirely in terms of the variables entering the many-body 
problem in the comoving frame. Hence we have transformed Eq. (I3.30p to a transparent 
physical form and demonstrated that it is indeed the correct equation of motion for the fluid 
elements. 

Using the results of Eqs. (I3.34p and (13.370 we can write down the following final system 
of equations describing the dynamics of the full many-body system 

idt^iii, ...,^N,t) = H[g^,, A]^{^u ...,^N,t) (3.38) 
-^itt) = ^^^\im{d^-d^>)pite,t) (3.39) 

mx^ = E^{^,t) + [i X B(x,t)]^ - (3.40) 

where the Hamiltonian H[g^y,A\, the reduced density matrix p[\E'](^, t), and the stress 
tensor P^y[gap, A,^'^] are defined after Eqs. (13.290 . (13.330 . and (I3.36p . respectively. The 
metric tensor g^u{^, t) entering Eqs. (I3.38P and (I3.40p is related to the Lagrangian trajectory 
via Eq. (Km . 

It is important to stress out that the system of Eqs. (I3.38l) - p.40p is mathematically 
equivalent to the original many-body Schrodinger equation (12. ip . In fact, everything we 
did to derive Eqs. (I3.38p - (l3.40p from Eq. (12. ip was an identical change of variables aimed 
at separating the convective and the "relative" motions of quantum particles. However 
after this identical transformation the structure of the many-body theory becomes quite 
remarkable. The key observation is that the physical external fields enter only the equation of 
motion for the fluid elements, Eq. (13.400 . while the many-body dynamics, which is governed 
by Eqs. (13.380 and (13.390 . depends only on the fundamental geometric characteristic of the 
comoving frame - the metric tensor 5'^i/(^,t). Formally Eqs. (I3.38p and (13.390 describe the 
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dynamics of quantum particles driven by a given time-dependent metric and constrained 
by the requirement of zero current density. This constrained many-body problem can be 
equivalently represented in a form of a Dirac-Frenkel variational principle with the folowing 
action functional 

S[^,A] = £dt (tmOtl^) - {^*\H[g^,,Am) . (3.41) 
Indeed, the conditions for the extremum of this action, 

are identical to Eqs. (13.381) and fl3.39p . respectively. From the variational formulation we 
clearly see the the effective vector potential is simply a Lagrange multiplier that ensures the 
zero current constraint. 

The formulation of the many-body part of the problem in a form of the variational princi- 
ple can be viewed as a time-dependent ^analog of the Levi-Lieb constrained search formulation 
of the static DFT 
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19|,[2^(see also [3l|). By finding an extremizer of the functional (13.411) . 
or, equivalently, by solving the constrained many-body problem of Eqs. (13.381) and (I3.39P for 
a given metric of the form (I2.25P we get the wave function \1/ and the Lagrange multiplier 
A as universal functionals of the metric tensor: \I' = \l'[5'^,y] and A = A[g^^i]. Substitution 
of these functionals in to Eq. (I3.36P gives the universal stress tensor functional P^ulQap]- 
Thus Eq. (I3.40p becomes a closed equation of motion for fluid elements, which determines 
the Lagrangian trajectories of the system. As all basic quantities are functionals of the met- 
ric tensor, which phy sically is nothing but the Green's deformation tensor of the classical 
elasticity theory 321], it is natural to call this approach the time-dependent deformation 
functional theory (TDDeFT). 

The direct solution of Eq. (I3.40p with a known functional P^u[gai3\{.$,ii) gives the de- 
scription of the convective motion in terms of the Lagrangian picture. Alternatively we can 
transform P^ulQapli^.i) to the laboratory frame to recover the tensor P^,^ [v] (x, t) , 

P^4v](x,t) = g!|i!i5^^[^„^(|(x,t),t)](|(x,t),t), (3.43) 

which can be used either in the hydrodynamic formulation of Eqs. (I2.7p - (l2.8p or to calculate 
the xc potentials for the KS formulation of TDCDFT, Eqs. (I2.17p - (l2.18p . Finally, since in 
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the laboratory frame the function ^(x, t) (the initial point of the trajectory that arrives to 
X at the time t) can be found from the equation 

[dt + v(x, t)9x]|(x, t) = 0, ^(x, 0) = X (3.44) 

the stress tensor determine! by Eq. fl3.43p is indeed a universal functional of the Eulerian 
velocity v(x, t). 

Therefore we recovered the full formal structure of the traditional TDCDFT, but at the 
fundamentally new level of understanding. First of all, in TDDefFT formalism a closed 
theory of convective motion appears from a regular and conceptually clean procedure: it 
is simply a first universal step in solving the many-body problem in the comoving frame. 
Hence now we clearly understand where the universal functional entering the theory come 
from and why they are universal. Also, the causal structure of all functionals becomes 
absolutely transparent. Indeed, both the constraint of Eq. (13.391) and the right hand side in 
Eq. (13.381) are local in time. Therefore the solution of the constrained problem of Eqs. (I3.38P 
and (I3.39p . provided it exists, defines the wave function as a strictly retarded functional 
of the metric tensor. Hence both P^ulgais] and -P^zy[v] depend, respectively, on ga/s and gag 
in a strictly retarded manner. Thus the new formulation completely removes all doubts 
concerning the predictable power of both hydrodynamic and KS formulation of TDCDFT. 



C. Basic uniqueness and existence theorems of TDDefFT 

The traditional formulation of TDDFT/TDCDFT is based on two key mathematical 
statements: (i) the uniqueness and invertibility of the mapping from external potentials to 
the den.ty/cu.e„t, and (.„ the ™jpte.entabUity of the den.ty/cuttent. At the pte^ent 



v-repres 

m\ 



stage of the theory we have proofs |3|, U, ll5|] of the mapping theorem for Taylor expandable 
potentials, while the v-representability problem, strictly speaking, remains open. To prove 
the v-representability within techniques used in the available proofs of the mapping theorems 
on has to demonstrate that the uniquely constructed power series for potentials do converge, 
which has not been done up to now. 

Within the new formulation based on TDDefFT there is no question of mapping as there 
is no external potential in the universal many-body problem in the comoving frame. Instead 
the two above mathematical issues reappear in a form of the uniqueness and existence 
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of a solution to nonlinear system of Eqs. (13.381) . (13.391) . The uniqueness is a reminiscent 
of the potential-to-current mapping, while the problem of the existence of a solution is 
equivalent to the problem of interacting v-representability. Equations (13.381) and (I3.39P 
correspond to interacting system, but the same set of questions can be asked for a system of 
noninteracting particles. In the later case the dynamics is characterized by one particle 
orbitals (fj{^,t), j = 1 . . .N, and the universal constrained many-body problem simplifies 
as follows 

idtipj = g'^id^^^ + A,,) {id^^ + A^)g~^ipj, (3.45) 

. N 

= ^ Y.^V*AVj - VAV]) (3.46) 

where no(^) = lv^j(^' 0)P — Iv^il^? density of particles which by 

construction is independent of time. The problem of existence of a solution to the sys- 
tem of Eqs. (I3.45p . (13.461) is similar to the well known problem of the noninteracting v- 
represent ability in the traditional formulation of the TDCDFT. Note that for the practical 
purpose of TDDefFT/TDCDFT it is not necessarily to prove the well-posedness of the sys- 
tems Eqs. (13.381) . (I3.39P and Eqs. (I3.45p . (13.460 for any metric. It it enough to consider 
the metric tensors of the form (I2.25p . i. e. the metrics which are generated by an invertible 
transformation coordinates and correspond to a flat space. 

The key advantage the new formulation of old issues is that now the basic questions 
underlying TDCDFT are posed in the standard form common in mathematical physics. 
It is quite likely that the uniqueness and the existence theorems for nonlinear systems of 
Eqs. (13.381) . (13.390 and Eqs. (13.451) . (13.460 can be proved using standard methods of the 
functional analysis and the theory of differential equations (see, e. g., 33|] and Chapter X in 

uUy applied to the analysis of the time-dependent Hartree 



Ref. |3j]), which have been success: 
and Hartree-Fock systems 
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37|. Structurally the nonlinear systems of Eqs. (13.380 . 



(I3.39P and Eqs. (I3.45p . (I3.46P are of the type of a time-dependent Schrodinger equation with 
a special cubic nonlinearity, which, to the best of my knowledge, have not been considered 
before. An encouraging observation is that, in contrast to most known systems nonlinear 
Schrodinger equations, our equations are exactly integrable in the simple case of A^ = 1 
(when the systems of Eqs. f l338|) . fl339|) and Eqs. ([335]), fl3:i6|l coinside) (l7|. The exphcit 
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solution solution for = 1 is 



^{tt) = VM^e'''^^^'\ A,{$,t) = d^MU) (3.47) 

where 

1 /"* r 1 

^ = MO + -r, ^ / 9'-'d^,^g^''d^.g--^^, dt' (3.48) 

and foiO is a phase of the initial state. The corresponding stress tensor functional, 
Eq. fl3.36p . takes the form 

1 r_ 1 , „ 1 , ^^aPf, ^0 



W = - {di.g-^^^m'^g--^^) - ^d^.^gg-^d^,^\ (3.49) 



Equations f l3.47p and fl3.48p clearly show that the wave function and the effective vector 
potential are retarded functionals of the metric tensor. The stress tensor for the one particle 
case turns out to be absolutely local in time. Adding more particles should make the 
time dependence nonlocal (though necessarily retarded), but hopefully it will not spoil the 
existence and uniqueness of a solution. 

Apparently, a Taylor expansion-based proof of uniqueness also works in the new setting, 
and actually becomes almost trivial. To show this we first substitute the constraint fl3.39p 
into Eq. fl3.38p to get a closed nonlinear evolution equation of the form 

%dt^ = E[g^,,^^ (3.50) 

with a local in time nonlinear Hamiltonian. Let us assume that the metric tensor g^^ is 
v-representable. In other words, it is given by Eq. fl2.25p with a Lagrangian trajectory that 
corresponds to a convective motion in some external potential. This is equivalent to the 



assumption that the solution to Eq. fl3.50p does exists. If in addition we assume that the 
metric possesses a covergent Taylor expansion in time, we can represent both gfj,u{t) and 
\t'(t) in a form of power series 



k=l k=l 



and insert them into Eq. f l3.50p . Since the right hand side of Eq. f l3.50p is local in time all 
the coefficients \E'*-'^-' with any k are trivially expressible by recursion in terms of \E'o and g^il 
with / < k. The later property is another direct manifestation of the retarded character of 
the functional \E'[(7^i,](t). Thus we have proved that for a given v-representable metric the 
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solution is unique and has correct causal properties required for TDDefFT/TDCDFT to be 
a predictable theory in any formulation. 

If we relax the v-representability assumption we would have to prove that the uniquely 
constructed power series for ^(t) is convergent, which seems to be technically a very hard 
task. Most likely, early or later, the existence theorem will be proved by more advanced 
mathematical methods, applications of which becomes possible within the present reformu- 
lation of the theory. It is also quite likely that the physically unmotivated requirement of 
the Taylor expandability will be removed, as it can be done in the linear responce regime of 
TDDFT [381. 



IV. CONCLUSION 



In this work we proved the equivalence of two approaches to the convective dynamics of a 
general quantum many-body system and reconstructed TDCDFT on the new grounds. The 
most important outcome is the possibility to formulate TDCDFT within a conceptually clean 
and straightforward two step procedure which resembles the constrained search formulation 
of the static DFT. On the first step one solves a constrained time-dependent many-body 
problem to find a stress tensor as a universal functional of the deformation tensor. On the 
second step we use this functional to calculate the evolution of the current and the density for 
a given configuration of external fields. In this formulation the the vector potential-current 
density mapping theorem, and the v-representability problem are restated as the uniqueness 
and the existence theorems for a solution of a certain time-dependent nonlinear Schrodinger 
equation. Hence the fundamental questions of TDDFT appear in a standard setting of 
mathematical physics. To my knowledge a very special type of nonlinearity appearing in 
TDDefFT has not been studied before. One of the hopes is that the new restatement of 
the basic theorems combined with a growing practical popularity of TDDFT/TDCDFT will 
attract attention of mathematicians and mathematical physicists to the formal foundations 
of the theory. I also believe that a clear physical picture behind the present approach should 
stimulate further developments of the theory within more physically oriented part of the 
community. In any case TDDFT and TDCDFT is still an exciting, active and promising 
area of research. Clearly, there is still much to be done both to put the theory on really 
firm mathematical grounds, and to develop practically working approximate functionalfor 
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numerous applications. 
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